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Abstract — The use of Ecosystem and Biodiversity 
mapping, land use land cover change detection has been 
advocated in preparation of developmental master plan in 
towns and cities. Noticeable changes have been observed 
within Malete Elemere community since the establishment 
of Kwara State University Malete, yet its spatial pattern 
and socio ecological implication have not been 
investigated. This work seek to determine and produce 
land cover land use change map of Malete Elemere over 
the last 10 years and post 15 year periods through change 
detection techniques so as to evaluate the impact of the 
establishment of Kwara State university on the settlement 
spatial development. Landsat 7 Enhanced Thematic 
Mapper Plus (ETM+) satellite images of 2005, 2010 and 
2015 of the study area were acquired from USGS at 
spatial resolution of 30 m. Radiometric correction were 
applied to all the images using radiance modules in 
Idrisi32 with radiance spectral value set at DN 0 (Lmin) 
and 255 (Lmax). An unsupervised classification was 
carried out on the composite images of bands 4,3,2,1 for 
all the selected years to identify possible maximum 
spectral reflectance classes, this was followed by 
supervised classification using training sample from the 
field survey from which image to image spatio-temporal 
changes statistics were extracted. To generate a 
prediction of LULC changes for 2025, Cellular 
Automata-Markovian transition estimator (CA-Markov) 
in Idrisi32 was used. Various Kappa statistics was used to 
evaluate the performance of prediction with an average K 
statistics of above 0.83 recorded. The result shows that 
built up area gained an astronomical increase (180%) 
between 2005 and 2015 while forest lost significantly 
(34%) within the same periods, with most of the gains 
occurring in 2010 and 2015 after the establishment of 


KWASU. By 2025, two Major growth pole centres will 
emerge along Malete Elemere Axis and one minor in 
Jenkunu Omoni Axis which will exert a great stress on 
infrastructural facilities and may create a chaotic 
condition if left unattended to. 

Keywords — Land use land cover (LULC) change, 
Spatio temporal, prediction, developmental planning. 

I. INTRODUCTION AND LITERATURE 
REVIEW 

Ecosystem and Biodiversity Mapping (EBM) has been a 
veritable tool being used by environmental managers and 
scientists for sustainable land use and planning of natural 
resources (Fuller et al 2014, Barthlother al. , 1999). EBM 
does not only provide information on spatial distribution 
of species across the landscape but also serve as vital 
source of information on species natural habitat, species 
values and functions, the level and magnitude of any 
disturbance in the ecosystem (land cover land use change) 
all of which have great implication on developmental 
planning (Hegazy and Kaloop, 2015). Given the rate of 
deforestation and loss of biodiversity especially in 
developing countries through carelessness, poor planning 
and high level of poverty which has put undue pressure on 
natural resources, it is practically challenging to attain 
sustainable development without adequate information on 
the ecosystem and the biotic and abiotic composition 
(Gladstone and Thomas, 1990). The use of ecosystem 
land use land cover change detection and biodiversity 
mapping have been advocated in preparation of 
developmental master plan in towns and cities. This could 
help development planners in identifying protected areas, 
open space and designing of zoning (BRC, 2013). The 
Biodiversity Resources Centre, New York United States 
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had undertaken the habitat (biodiversity) mapping project 
for over ten towns in Hudsonia developed areas as a tool 
for town and country planning. 

Geographical Information System and Remote Sensing 
(GIS/RS) have proven to be very useful for large scale 
mapping of ecosystem and land cover (Trisurat^ al ., 
2000; Foody 2002; Lu and Weng, 2007). These 
approaches are faster and enable wider geographic 
coverage within limited time frame (USGS/GAP, 2002; 
Lowry et al., 2005). Many studies on land cover and 
vegetation/ecosystem mapping have used data from 
Advanced Very High Resolution Radiometer AVHRR 
(Defries and Townsend, 2002), Multispectral Moderate 
Resolution Imaging Spectra Radiometer-MODIS (Xiao et 
al., 2002), and Landsat Enhanced Thematic Mapper Plus- 
ETM+ (Lu and Weng, 2007; Yuan et al., 2006; Yang and 
Lo, 2002). While, AVHRR was originally designed for 
meteorological service and has only two spectral bands- 
red and near infrared which although sufficient for basic 
vegetation study, MODIS though has low spatial 
resolution yet has more spectral bands including short 
wave infrared (SWIR) which can be used for obtaining 
greater details and advanced vegetation analysis such as 
leaf moisture, soil moisture, canopy water contents among 
others (Boles et al., 2006;Cacceto et al., 2002a, b) 
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Malete and its adjoining settlements were and still are 
rural communities with the establishment of Kwara State 
University Malete (KWASU) campus in 2009. It has 
since witnessed significant physical infrastructural 
development many of which are done with little or no 
consideration for its ecological implication, now that the 
development is still at its early stage and given the vision 
of KWASU to prepare a Development Master Plan. 

1.1 Aim 

To assess the landscape dynamics prior and since the 
establishment of KWASU and predict the socio economic 
and ecological implication on the adjoining community. 

1.1.1 Objectives: 

- To assess the landscape dynamics prior to and since 
the establishment of KWASU and predict the socio 
economic and ecological implication on the 
adjoining community. 

- To determine biodiversity loss/gain over 15 year 
periods through change detection techniques and 
highlight its implication on developmental Planning. 

- To predict possible land use pattern in the next 15 
years and the relevant planning strategies to adopt. 

II. MATERIALS AND METHOD 

2.1 Study Area 


4'IOtTE 4-20XTE 4'30'CTE 4°400’E 4’50'0"E 



•0"N- 


■9°0'0 


4°25'30’E 4”27'30"E 


Fig.l: The study area with Nigeria map inset 
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The study area is located in Moro Local Government 
Area of Kwara State and lies within latitude 8.6563°N to 
8.8136°N and longitudes 4.2359°E to 4.5410°E. It 
comprises of Malete, Elemere, KWASU Campus, and the 
adjoining communities covering an area of about 
(157,701 Hectares) of land. 
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The study area is about 25 km North of Ilorin, the Kwara 
State capital though a relatively virgin area, it is highly 
vulnerable to unplanned expansions due to its proximity 
to the state capital and recently the siting of KWASU 
campus. 

2.2 Methodology Flow Chart 





r 



l 



Fig. 2: Methodology flow chart 


2.3 Materials 

The geographic extent of the study area was first 
determined and the shapefile prepared in Arc GIS 10.2. 
All the satellite images were co registered to the same 
study area shapefile to give similar spatial dimension. 


Archive Landsat 7 Enhanced Thematic Mapper Pus 
(ETM+) satellite images of 2005, 2010 and 2015 of the 
study area were acquired from USGS at spatial resolution 
of 30m. 


Table. 1: Summary of data acquired and used 


S/N 

Data Type 

Years of 
Acquisition 

Resolution / 

Scale 

Source 

Application 

1 

Landsat ETM+ 

Images 

2005, 2010 

2015 

30m 

GLCF/USGS 

Land cover Mapping 
(NDVI) 

2 

Digital Elevation 

Model (DEM) 

2000 

90m 

SRTM/NASA 

Physiographic Mapping 

3 

GeoEye-1 

2016 

40cm 

TerraServer 

Image 

V alidation/Classification 

4 

Climatic Data 

1900-2000 

2000-2013 

0.86km 2 

WorldClim / 

URBDA 

GPS/GBIF 

Bioclimatic classification 

5 

Species Occurrence 

2005 

Nil 

Google Earth 
Field work 

Suitability index 

6 

Soil Data 

2007 

250m 

UNEP 

Physiographic Mapping 

7 

Vegetation data 

2013 

25mx4m 

Field Work 

Veg/Classification 

Acronyms: 

Enhanced Thematic Mapper Plus 

= ETM+ 
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Global land Cover facility = GLCF 

Shuttle Radar Topographic Mission = SRTM 

Global Biodiversity Information Facility = GBIF 

Upper Benue River Basin Development Authority = URBDA 

National Aeronautics and Space Administration = NASA 

United Nation Environmental Programme = UNEP 


2.4 Image Pre-Processing 

Radiometric corrections were applied to all the images 
using radiance modules in Idrisi 32 with radiance spectral 
value set at DN 0 (Lmin) and 255 (Lmax). All images 
were collected in the months of July respectively being 
the possible periods of rainy seasons for effective 
measurement of plant vigour. Attempt was made to 
collect cloud free images in all the time series. Images 
band combinations were performed on bands 3, 2, 1 and 
4, 3, 2 for classification (urban, water bodies and 
agriculture) and vegetation differencing (forest and 
grassland) respectively. 

2.5 Image Processing 

The unsupervised classification was carried out on the 
composite images of bands 4,3,2,1 based on pixel spectral 
characteristics/signatures of various land cover. An 
iterative ISODATA (Maximum Likelihood Classifier) 
algorithm was used in Arc View GIS (version 10.2) and 
Multispec (2013 version). This grouped similar pixels in 
the image into clusters or categories, and help us in 
determining maximum spectral classes in the images. 
There was no significant class change after eight spectral 
classes thus provided a good idea on possible classes for 
our classification. 

To enhance our classification and identify the green index 
or plant cover in the study area, the Normalized 
Difference Vegetation Index (NDVI) was performed to 
compliment the earlier unsupervised classification. NDVI 
is a remote sensing /GIS techniques used over the years 
by scientists to quantitatively and qualitatively evaluate 
the vegetation covers of an area (Neelima T.L et al 2013). 
NDVI as proposed by Rouse, et al (1974) is 
mathematically defined as: 

NIR - R 


Where, NIR and R are the reflectance in the near infrared 
and red regions respectively. It is the algebraic 
combination of red and near infrared bands to represent 
the amount of green vegetation in the image. In the 
NDVI, the values for a given pixel value is always in a 
number that ranges from -1 to +1. A zero means no 
vegetation and close to 1 indicates the highest possibility 
of green leaves (Biehl, 2010). 

Field survey/ground truthing of 120 points were 
conducted on 20-30th July 2016 and was used as our 


training sites. This were overlaid on unsupervised land 
cover classes and combined with GeoEye-1 images of 
2016 to prepare the supervised classification (Salako et 
al., 2016). The result of supervised classification 
produced the following classes of land use land cover in 
the study area: Forest, Mixed forest, Grassland, Farmland, 
Adjoining built up and Built up. 

2.6 Land use land covers projection: CA Markov 
Techniques 

To generate a prediction LULC changes for 2025, the 
Markovian transition estimator in Idrisi32 was used. 2005 
land cover image was input as earlier image while 2015 
land cover image was used as later or second image with 
the number of time periods between the first (earlier 
image ) and second (later image ) was 10 while the 
number of time for projection from the second image was 
also set at 10 years that is 2025. Equal probability was 
assigned to the entire pixel under estimation. Based on 
this the following estimation was generated: the 
probability transition matrix (Table 2), the transition area 
matrix (Table 3) and conditional probability image. To 
add the spatial dimension to our prediction the cellular 
automation (CA) was combined with Markov transition 
estimation with 2015 land cover image used as basis for 
projection and the earlier generated transition area matrix. 
The cellular automation was set at 10 to project for 2025. 

2.7 Model Validation 

An important stage in the development of any predictive 
change model is validation. Typically, one gauges one 
understanding of the process, and the power of the model 
by using it to predict some period of time when the 
landcover conditions are known. This is then used as a 
test for validation. IDRISI supplies a pair of modules to 
assist in the validation process. The first is called 
VALIDATE, and provides a comparative analysis on the 
basis of the Kappa Index of Agreement. Kappa is 
essentially a statement of proportional accuracy, adjusted 
for chance agreement. However, unlike the traditional 
Kappa statistic, VALIDATE breaks the validation down 
into several components, each with a special form of 
Kappa or associated statistic based on the work of Pontius 
(2000): Kappa for no information = Kno • Kappa for 
location = Klocation • Kappa for quantity = Kquantity • 
Kappa standard = Kstandard • Value of Perfect 
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Information of Location = VPIL • Value of Perfect 
Information of Quantity = VPIQ With such a breakdown, 
for example, it is possible to assess the success with 
which one is able to specify the location of change versus 
the quantity of change. The accuracy of prediction is 
measured by the performance of various K statistics, the 
higher the value the better the prediction. 

III. RESULT AND DISCUSSION 
3.1 LULC 2005 - 2015 

In 2005 forest cover (mixed and closed) constituted over 
32 % of total LULC with built up covering only 876 ha of 
land representing 6% (Fig.3). While, open but adjoining 
built up constituted about 11%. Five years later, in 2010, 
part of adjoining built up had changed to cropland thus 
cropland increased from 19% in 2005 to 25% due to 
Fadama project established at the period (Fig. 4). The 
Built up was almost at stable point remaining at 5%. By 
2015 the effect of siting KWASU campus in Malete has 
become obvious, the built up had increased from 5% in 
2010 to 15% in 2015 covering 2400.03 ha. of land from 
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803 ha in 2010 (fig.5). Many estate developers and 
private builders sought for land within KWASU campus 
and Malete . However forest cover has been worst hit It 
fell from 17% in 2005 to about 11% in 2015. This was 
noticeable in Western Bi Ala where about 150 ha of forest 
land changed to shrubby forest and grassland in 2010 and 
by 2015 reduced to a narrow strip of forest of less than 45 
ha (fig. 5) 

Percentage of change analysis between 2005 and 2015 
revealed that adjoining lowland and forest cover were the 
top losers with about 66% of adjoining lowland lost to 
either cropland and or built up (Fig. 6). This was followed 
by forest cover which lost about 34% of their total land 
(1500 ha). The top gainer was the built up area which 
recorded an astronomical increase of 180 % totalling over 
1600 hectares of land (Table 2). This was noticeable at 
the major settlements of Malete, Elemere, and KWASU 
campus with several residential buildings used either for 
student hostels or private residences within 500 m radius 
of the campus. Smaller settlements like Apodu, Jenkunu 
and Gbugudu increased by 25% between 2009 and 2015. 


Table.2: LULC change analysis 2005- 2015 


Classes 


Hectares 

2005 


% 


Hectares 

2010 


% 


Hectares 

2015 


% 


Open /Built up 

876.42 

5.5 

803.25 

5.2 

2400.03 

15.4 

Adjoining built up 

1695.96 

10.9 

105.57 

0.7 

568.17 

3.7 

Cropland 

2986.83 

19.2 

3796.74 

24.6 

4529.16 

29.1 

Grassland 

4879.71 

31.4 

3782.43 

24.5 

3210.12 

20.6 

Mixed forest 

2529.09 

16.3 

4862.97 

31.5 

3142.08 

20.2 

Closed forest 

2591.73 

16.7 

2075.22 

13.5 

1710.18 

11 

Total 

15559.74 

100 

15426.18 

100 

15559.74 

100 



Legend 

• Major settlements 
HH Open/built up 

] Built up adjoining land 
| Cropland 


Grassland 
Mixed forest 
Closed forest 



Fig. 3: LULC 2005 
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Fig.4: LULC2010 



Fig.5: LULC 2015 
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Fig.6: LULC 2005- 2015chart 


Efllt up, 180 

150 200 


www.ijeab.com 


Page I 2093 


















International Journal of Environment, Agriculture and Biotechnology (IJEAB) 
http ://dx.doi. org/10.22161/ijeab/2.4.6 7 


3.2 Probability transition matrix for 2025 

Probability matrix of land use land cover change for 2025 
was generated using CA Markov model (Table 3). Most 
land use tend to transit to cropland and Mixed forest 
especially in Elemere segment, the two LU had the 
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highest transition probability matrix of 0.25 as compared 
to the overall average of 0.13 The predicted land use 
changes are as follows: Open adjacent to Built up, Forest 
to Mixed Forest, Grassland to Cropland. 


Table. 3: Probability matrix 2015 and 2025 

2025 



OPB 

ADJ 

CRP 

GRS 

MF 

FR 

OPB 

0.0192 

0.0871 

0.2082 

0.1733 

0.2853 

0.2269 

ADJ 

0.0312 

0.1136 

0.2331 

0.1836 

0.2547 

0.1838 

CRP 

0.0473 

0.1449 

0.2556 

0.1911 

0.2197 

0.1413 

GRS 

0.0646 

0.1737 

0.2703 

0.1941 

0.1887 

0.1086 

MF 

0.0859 

0.2022 

0.2760 

0.1922 

0.1597 

0.0840 

FR 

0.1250 

0.2416 

0.2702 

0.1836 

0.1216 

0.0581 


OPB= Open and Built up area GRS= Grassland 

ADJ= Land adjacent to Built up MF= Shrubby/mixed forest 

CRP= Crop/farmland FR= Closed/dense forest 


3.3 Land Use Land Cover Projection 2025 

Land use Land cover change for 2025 was done using CA 
Markov model. This explains the probability transition 
matrix and area change calculated from ArcGIS and 
Multispec. The pattern observed between 2005 and 2015 
persisted in the projection with built up predicted to be 
having the higher percentage gain in land cover land use 


statistics of about 32% by closing up the adjoining open 
land (Fig. 8) while crop land especially around Elemere 
had a gain of over 4 % Mixed forest equally rose close to 
5% by 2025 and this was probably due to forest 
degradation which would be losing over 17% of its 2015 
hectares of land (1710.2) to about 1408.14 in 2025 
(Table 4). 



Fig. 7: Projected LULC 2025 


Classes 


Built up 

Cropland/Farmland 
Grassland 
Shrubby/Mixed Fore 
Closed Forest 


Table.4: Percentage change analysis of LULC 2015 2025 


Class 

2025 


2015 


2015/2025 


Area (Ha) 

% 

Area 

% 

A % 

Forest 

1408.14 

9.1 

1710.2 

ii 

-17.2 

Mixed Forest 

3268.71 

21.2 

3142.4 

20.2 

4.95 

Grassland 

2572.02 

16.7 

3210.1 

20.6 

-18.91 
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Crop/Farmland 

4683.63 

30.4 

4528.3 

29.1 

4.46 

Adj Built up 

520.69 

2.3 

567.5 

3.7 

-37.83 

Built up 

3066.11 

20.3 

2400.1 

15.4 

31.82 

Total 

15519.3 

100 

15558.6 

100 



2015/2025 %A 

■ Forest Mixed Forest Grassland ■ Farmland Adj Built up ■ Built up 




Built up 




Adj Built up 


Farmland 

4.46 

-18.91 Grassland 




Mixed Forest 

4.95 



-17.2 Forest 



Fig.8: LULC Changes 2015 and 2025 


3.4 Model Evaluation 

For validation of Markov CA model using various Kappa 
statistics, 2015 observed image was used as referenced 
and was compared with the simulated image of 2015 to 
see the similarity between the actual and projected land 
use land cover map the following Kappa statistics was 


generated (Table 4). The result shows high performance 
of the model and its prediction for 2025 and 2030 since 
most K Statistics (K standatd value of about 0.893 and 
Klocation of 0.922 Kno -0.8937) were above 80% 
(Praveen-Subedier al. , 2013). 


Table.5: Validation of projected 2015 LULC map with actual 2015 LUC map 




Ability to Specify Quantity 



No[n] 

Medium[m] 

Perfect[p] 


Perfect[P(x)] 

P(n) = 0.9007 

P(m) = 0.9754 

P(p)= 1.0000 

Ability to Specify 
Location 

Medium[M(x)] 

M(n) = 0.8461 

M(m) = 0.9150 

M(p) = 0.9381 


No[N(x)] 

N(n) = 0.2000 

N(m) = 0.1996 

N(p) = 0.2053 

CorrectChance = 0.1996 

CorrectQuantity = 0.0000 

CorrectLocation = 0.7154 

ErrorLocation = 0.0604 


ErrorQuantity = 0.0246 PerfectChance = 0.2000 PerfectLocation = 0.7947 PerfectQuantity = 0.0053 VPIL = 
0.0604 VPIQ = 0.0231 Kno = 0.8937 Klocation = 0.9221 Kquantity = 0.7488 Kstandard = 0.8938 


IV. CONCLUSION AND RECOMMENDATION 

• The adjoining built up will be closed up and fully 
merged with the built up area for developmental 
project by 2030. 

• Malete-Elemere growth pole axis is at risk of 
chaotic urban growth if action is not taken now. 


BialaBudo Are forest will be lost and transit to 
shrubby forest and perhaps grassland thus deplete 
the area of high biodiversity values. 

The North East section is potentially forest reserve 
zone and could be designated as conservation area. 
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• Sustainable land use plan should give high priority 
to enhancement of community/land biodiversity 
value. 

• Full community participation and input in land use 
planning. 

• High level of adaptability in response to change 
with emphasis on bottom top approach to 
planning. 


• Planning based on up-to-date data and full 
integration of geospatial data. 

• Functional government agency to regulate, 
administer and implement plan policy e.g. Malete 
Elemere Development Area Commission. 

Recommendation: Three zones of different 
land use planning are recommended 


3 Zones different needs 



Elemere=Agric. semi urban 


ClHMI 

■ r«t«t 

Ubortfwrtf 

GnttM 

CnffiHwimd 



Malete Forest Agric. Urban 
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